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ANALYSIS OF A LATERAL PI LOT- INDUCED 


OSCILLATION EXPERIENCED ON THE 
FIRST FLIGHT OF THE YF-16 AIRCRAFT 
John W. Smith 

Dryden Flight Research Center 


INTRODUCTION 


During a first flight, or early in the flight test program of new aircraft, 
numerous incidents of pilot induced oscillations (PIO's) have occurred (ref. 1). 
In a majority of the cases, the occurrence has been about the longitudinal axis, 
such as in the experiences reported and analyzed in reference 2. However, the 
pilot can also develop and sustain a PIO in maintaining roll attitude control. 

During a high speed taxi run, scheduled as part of the buildup prior to 
the first flight, the YF-16 aircraft inadvertently became completely airborne; 
this necessitated a complete go-around. The test plan called for a taxi run to 
peak speed, and a minimum lift-off, to be followed by a slight lateral control 
input. However, a fairly high frequency PIO (.7 cycle per second) occurred on 
lift-offand was sustained for 12 complete cycles. As a result, the aircraft 
became misalined with the remaining runway. The pilot, therefore, was compel- 
led to add power and complete the takeoff (appendix A). The magnitude and rate 
of the pilot roll control inputs were sufficient to position and rate saturate 
the roll control system. 

The PIO occurrence was reported by the Air Force in reference 3. This 
experience, plus the studies conducted by the contractor, led to a system 
design change which in general reduced the overall lateral control system gain 
for the takeoff and landing phases. The interim system that was implemented, 
following this first flight, is identified in this report as the modified flight 
control system (MFCS). 

Despite the long history of PIO's, a generally valid prediction technique 
does not exist. One of the means available to designers for their prevention 
is to compare the characteristics of new aircraft with the characteristics of 
aircraft that have proved to be resistant or prone to these oscillations. 

Because the YF-16 has a highly advanced fly-by-wire control system, a 
study of its PIO situation would contribute to the body of knowledge on air- 
craft PIO characteristics. Consequently a study was initiated to investigate 
the integrated control characteristics of the aircraft, control system, and 
pilot. The problem was approached and analyzed initially in a linear fashion 
with and without ground effects included. Nonlinear analytical studies were 



utilized to predict the frequency responses under saturated conditions and used 
to compare the initial control system to the MFCS control system with and with- 
out ground effects. 

This report discusses the theoretical analyses made at the Dryden 
Flight Research Center regarding this PIO occurrence on the YF-16 prototype air' 
plane. 


SYMBOLS AND ABBREVIATIONS 


Physical quantities in this report are given in the International System 
of Units (sp and parenthetically in U.S. Customary Units. The measurements 
were taken in U.S. Customary Units. Factors relating the two systems are 
presented in reference 4. 

A amplitude 

A/C aircraft 

A/C(s) aircraft transfer function 

Aj^ actuator 


b 

C 


lateral acceleration, g 


reference span, m (ft) 
generalized control output 
rolling-moment coefficient, 


Rolling moment 
qsb 


C 


roll -damping derivative. 



rad 


-1 


C 


Z 


r 


rolling moment due to yaw rate. 




effective dihedral derivative, 



ap 


, deg 


-1 


2 



^9 -1 

horizontal tail-roll effectiveness derivative, deg 


flaperon-roll effectiveness derivative, v ■ , deg" 


9 -1 

rolling moment due to rudder deflection, deg 


^^9 -1 

rolling moment due to ground effect, rad 

<}> 


yawing-moment coefficient, 


Yawing moment 


-1 

yawing-moment due to roll rate, — rad'"^ 

9 2V 


yaw-damping derivative, — rad 

9 ^ 


-1 

directional -stability derivative, deg"'^ 


yawing moment due to aileron deflection, deg 


yawing moment due to flaperon deflection, 5 —^, deg “ 
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I u 


8C . 

rudder-effectiveness derivative, — 2^ deg”'^ 

86„ 


ac 

yawing moment due to ground effect, rad”"^ 


side- force coefficient, 

qs 


ac. 


V -1 

side-force derivative, deg 


aCw _i 

side-force derivative, deg 


SY -1 

side-force derivative, deg 

dp 


-1 

side force due to aileron deflection, deg 


3Cy _i 

side force due to flaperon deflection, deg ^ 


aCy 

side force due to rudder deflection, deg 


. aCy _i 

side-force derivative due to ground effect, deg 


mean aerodynamic chord, m (ft) 
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D 

DN 

e 

F 

F. 




f( ) 
G. 




g 

H(s) 

HT 

h 


'x> Iz 


‘xz 


denominator 

down 

error signal, rad or deg 
flaperon 

lateral stick force, kg (lb) 

pilot gain, kg/ rad (Ib/rad) 

function of parameter inside the parentheses 
horizontal tail forward loop transfer function 

flaperon forward loop transfer function 
rudder forward loop transfer function 
acceleration due to gravity, 9.8 m/sec^ (32.2 ft/sec^) 
roll rate feedback transfer function 

yaw rate feedback transfer function 

lateral acceleration feedback transfer function 

feedback transfer function 
horizontal tail 
altitude, m (ft) 

vehicle moments of inertia about the x- and z-body axes, 
respectively, kg-m^ (slug-ft^) 

product of inertias, kg-m^ (slug-ft^) 


jw 


imaginary part of Laplace transform variable, rad/sec 
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K 

KN 


LEF 

M 

MAC 

MFCS 

N 

PIO 

P 

q 

R 

r 

S 

SAS 

s 

T(s) 

TEF 

t 

V 

W 

X 

Y(s) 

a 


gain constant 

gain constant with a limit 

roll rate gain 

roll attitude gain 

leading edge flaps 

Mach number 

mean aerodynamic chord 

modified flight control system 

numerator 

pi lot- induced oscillation 

rolling angular rate, rad/sec 

dynamic pressure, N/m^ (Ib/ft^) 

generalized input 

yawing angular rate, rad/sec 

reference planform area, m^ (ft^) 

stability augmentation system 

Laplace transform variable, (a + jiu), rad/sec 

transfer function 

trailing edge flap 

time, sec 

true airspeed, m/sec (ft/sec) 
weight, kg (lb) 
input or output quantity 
generalized transfer function 
angle of attack, deg 
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p 

A 


angle of sideslip, deg 
transfer function denominator 


6, differential horizontal tail deflection, ^(6 mt/o " 6mt/i) 

^ deg or rad 

6p differential flaperon deflection, 

6^ rudder deflection, deg or rad 


t 

o 

<t» 

0/F, 

U) 

Subscripts: 

cr 

c 

e 

H 

HT 

L 


damping ratio of second-order response 

real part of Laplace transform variable, rad/sec 

angle of bank, rad 

system gain, rad/N (rad/lb) 

phase angle, deg 

frequency, rad/sec 

critical 
commanded 
error signal 
hydraul ic 
horizontal tail 
left 


n natural 

R right 

R generalized input 

X, Y, Z vehicle body axis 

1, 2, 3, sequential calculation in the order indicated. 

A dot over a quantity denotes the first derivative with respect to time. 
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VEHICLE DESCRIPTION 


Ai rcraft 

The YF“16 aircraft, under Air Force contract, was designed and manufac- 
tured by the General Dynamics Convair Aerospace Division of Fort Worth, Texas. 
The aircraft is a prototype lightweight fighter designed to have high g maneu- 
verability plus supersonic speed capability. Shown in figure 1 is a three-view 
drawing illustrating the relative size and location of the primary control sur- 
faces. The pertinent physical characteristics of the aircraft are presented in 
table 1. 


Primary and Secondary Flight Control Systems 

Aerodynamic surfaces . - Flightpath control in pitch and roll is achieved 
by a differential all-movable horizontal tail. Additional roll control is pro- 
vided by wing-mounted flaperons. Directional control is accomplished by a con- 
ventional rudder mounted on the trailing edge of a single vertical stabilizer. 

Fly by wire. - A full fly-by-wire system is utilized in all three axes. 

The system is quadruple redundant as to the sensors and the electrical output of 
the pilot controller, computers, and electronic circuitry. The mechanical 
inputs to control surface power actuators are provided by the electromechanical 
command servos and are single-fail-operate components. 

Secondary systems. - Leading-edge flaps are provided for takeoff, landing, 
and high angle of attack maneuvers. The leading-edge flaps are also scheduled 
as a function of angle of attack. Speed brakes located inboard of the horizon- 
tal tail provide deceleration capability. 

Pilot controls and augmentation. - Through a minimum displacement force 
sidearm controller, the pilot provides pitch and roll inputs to a rate command 
augmentation system, rather than controlling surface positions directly. Longi- 
tudinal force inputs command a blend of pitch rate and normal acceleration. 

The pitch feedback loops attempt to provide response characteristics that are 
almost invariant with flight condition. Lateral force commands roll rate propor- 
tional to the error signal (fig. 2). Favorable roll-yaw coordination is accom- 
plished through the angle of attack range by an aileron-rudder interconnect. 

A typical stability augmentation system (SAS) is implemented in the yaw 
axis and summed with rudder pedal position to provide both dynamic and static 
stability. Aircraft dynamic stability is augmented by a washed-out yaw rate 
feedback. Static directional stability is augmented by a lateral acceleration 
feedback. The surface authorities and rate limits are presented in table 2. 

A more complete description of the total flight control system and its 
design is given in reference 5. 
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Instrumentation 


The YF-16 instrumentation (ref. 3) was designed, installed, and maintained 
by G6nGrs1 Dynsniics* Th 0 ddts pr8S6nt6d in this rGport wers rGcordGd on bosrd 
by apulse code modulation system at 40 samples per second. The data system is 
considered to be accurate to within 3 percent of full-scale range of the param- 
eter being recorded. 


ANALYSIS AND DISCUSSION 


YF-16 PIO Experience 

Throughout the design and manufacturing phases of the prototype aircraft, 
the flight control systems underwent design, development, and testing on both 
fixed base and in-flight simulators. During these tests, the augmentation sys- 
tems were integrated with the basic aircraft response characteristics to pro- 
duce excellent flying qualities throughout the airplane's flight envelope. In 
addition, the feel characteristics (roll rate per stick force) selected and 
implemented were those felt to be the most desirable. However, as only natural, 
simulator studies motivate less attentive controlling or piloting than actually 
flying the aircraft. Therefore, it is highly likely that pilot gain can be 
much higher during flight testing than during any type of simulator flying 
especially on the first flight. 

Consequently, on the first flight of the YF-16 aircraft, which was sched- 
uled as a high speed taxi run, the pilot experienced a severe PIO immediately 
upon lift-off (appendix A). The magnitude and rate of the roll control input 
by the pilot was sufficient to position and rate saturate the roll control 
system. ^ The normal system lag and a constant lag due primarily to rate satu- 
ration introduced enough total system lag at the input frequency and amplitude 
to cause the pilot's lateral stick input/output relationship to be 180 degrees 
or more out of phase with the estimated peak bank angles. Table 3 lists the 
flight conditions of the aircraft during the PIO. 

Figure 3 presents a time history of the PIO experienced during takeoff on 
the first flight (from ref. 3). The parameters shown are lateral stick force, 
lateral control surface positions, and roll rate. As indicated by the satura- 
tion of the lateral stick force transducer trace, F , maximum roll rate was 

commanded. ^ The flaperons were the only control surfaces exhibiting continuous 
rate limiting during the PIO, as indicated by the triangular wave shape. The 
maximum rate for the sum of the two flaperon surfaces was approximately 112 
degrees per second (approximately twice the individual surface rate limits 
(tabled)), at a frequency of .67 cycle per second, and was sustained for 
approximately 12 complete cycles. For most of the cycles, the aircraft roll 
rates reached average peak values of about ±50 degrees per second, and bank 
angles were estimated to be ±15 degrees. Presented in table 4 is a summary of 
the pertinentcharacteristics of the PIO which were determined from an analysis 
of the time history. Also listed in the table are conditions common to other 
PIO's and unique features of the YF-16 system. 
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Roll Control Modeling 

The discussion that follows defines the modeling aspects and analytical 
description of the roll control mode for the YF-16 aircraft. 

System description and control elements. - Presented in figure 4 is a 
block diagram representation of the lateral control loops, including the pilot, 
that is suitable for both the linear and nonlinear PIO analysis. The pilot 
attempts to minimize the aircraft bank angle error to a wings- level position. 
The amount of force he applies per unit error signal is referred to as pilot 
gain. The original system roll command force gradient, dual lag filter, and 
actuator response descriptions are presented in figures 5 to 7. The roll stick 
force command gradient is described by a cubic equation shown by the dashed 
line in figure 5. The dual lag command prefilter consists of four nonlineari- 
ties (fig. 6), two in each feedback path. In each feedback path and in between 
the nonlinearities is a washout network which allows the affected loop to null 
when commands have a constant value. The purpose of this filter is to command 
a faster response initially and then a somewhat slower response as the error 
signal (e) approaches zero. However, the average phase lag per cycle can be 

closely matched by just the forward loop transfer function ( — if sine 

wave amplitudes limited to that of the saturated input are used in the mathema- 
tical model. 

At a flap setting of 20 degrees, the flaperons are at the full down 
deflection limit. A flap crossfeed logic utilizing the limiters and thresholds 
(fig. 2) maintains the gain of the roll control channel by doubling the deflec- 
tion of the upward moving flaperon to compensate for the lack of a downward- 
moving flaperon on the opposite side. Also, whenever the surfaces were moving 
in opposite directions, the total rate limiting would be twice the rate limit- 
ing exhibited by a single surface (approximately 112 deg/sec). A diagram typi- 
cal of the mechanical servos and actuators is presented in figure 7. The sys- 
tem rates, authorities, and servo characteristics for each system affecting 
roll response are presented in table 2. Rate limiting can be approximated with 
reference to figure 7 by using the following system. The restraint of any gain 
change is considered to be just in the forward loop of power actuator servos, 
and not in the feedback loop. 


R 
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or C, the output, would be 

C = N X I 

where the rate restriction N is 


Since the change in feedback gain due to the nonlinearity is neglected, this 
makes the closed-loop input/output relationship for saturated inputs slightly 
smaller than the actual system at the lower frequencies of interest (lu < 8 
rad/sec). 

The differential horizontal tail was unaffected, since surface limiting 
due to roll command in that path was not prevalent during the PIO. 

Modified flight control system . - Following the PIO experienced on the 
first flight, the YF-16 aircraft's roll control system was modified to include 
a takeoff and land mode. In general, these modifications reduce the system 
gain in the forward loop, reduce the roll command force gradient, and lower 
the maximum roll rate that can be commanded. In this report, this system is 
designated as the modified flight control system (MFCS). 

Presented in figure 8 is the roll command gradient implemented on the MFCS 
The dashed line indicates the cubic equation in the subsequent studies that was 
used to define the roll rate commanded per lateral stick force. Maximum com- 
manded roll rate was limited to 167 degrees per second. The forward loop gain, 
KN, of the roll control system was reduced from 0.20 for the original system to 
0.12 radian/radian per second for the MFCS (fig. 4). 

Aerodynamic characteristics. - The aerodynamic stability and control deri- 
vatives were obtained from rigid model tests conducted in wind tunnels at vari- 
ous facilities. The data were compiled and reported in reference 6. Table 5 
lists the pertinent coefficients as transferred to the body axis. The transfer 
functions ratioing roll rate to surface deflection were obtained by using typi- 
cal computer programs and the appropriate geometry and flight conditions 
(tables 1 and 3). 

Ground effects. - Airplanes flying near the ground plane during takeoff 
and landing experience changes in aerodynamics forces and moments that are 
termed ground effects. Many studies have been made regarding these effects, 
but they have addressed primarily the changes in forces and moments about the 
longitudinal axis. The studies have shown that near the ground plane the over- 
all effect is an increase in lift, a decrease in pitching moment, and a reduc- 
tion in drag force. 
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Lateral and directional changes in moments are for the most part con- 
sidered small, because the bank angle is usually near zero. However, if the 
bank angle becomes relatively large (cj) = 15 deg or greater), the change in 
pressure distribution may be large enough to affect the lateral dynamics. The 
extent of these effects, particularly regarding the flying qualities, would, 
of course, depend on the type and amount of stability augmentation. 

To complement these studies, and because this PIO was in the ground effect 
regime, an estimate was inserted in the equations of motion to account for these 
changes in the lateral forces and moments. 

Datcom (ref. 7), which is a compilation of design information previously 
published in theoretical reports, presents a section on methods of estimating 
the changes in longitudinal forces and moments as a function of height above 
the ground level. In this report, the estimates utilize the same information 
except that bank angle is related to the height difference between the right 
and left wing quarter-chord-of-the-MAC. Thus, the difference in lift force 
would cause a resultant rolling and yawing moment component. For a bank angle 
of 15 degrees and an average height of 2.74 meters (9 feet), it was estimated 
that: 


= -.026 per radian 
= .00635 per radian 
Cw = 1.376 per radian 

The transfer functions including- these additional terms are presented along 
with the basic aerodynamics in table 5. 

As a check, was calculated in a more straightforward manner with the 

the digital computer. The computer program consisted of a numerical computa- 
tion, utilizing the theory presented in reference 8, that resolved the forces 
on a lifting wing in the presence of a ground plane. The calculation for the 
same height predicted a smaller value = -.02 per radian). However, the 

larger value was used in deriving the transfer functions, to be conservative 
in assessing the influence of ground effect. 


Integrated Analysis of the PIO 

Linear calculation. - Pilot-in-the-loop, closed loop transfer functions 
were derived from aircraft transfer functions (listed in table 5) and the sys- 
tem gains. The closed- loop polynomial ratios were then used to predict the 
dynamic characteristics of the lateral flight control system for inputs below 
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position and rate saturation. Initial assessments of ground effect were inves" 
ti gated comparatively in this fashion. The linear calculations were also used 
as a special check case for subsequent nonlinear calculations. However, in 
deriving the pilot closed- loop response, some generalization is required regard- 
ing the feel system because of the nonlinearity associated with commanded roll 
rate and stick force. For convenience, the outer loop gain is considered to be 
the pilot gain times the feel gradient. 

Figure 9 makes a root locus (pilot closed-loop) comparison of the initial 
control system implemented on the prototype aircraft in and out of ground 
effects. The transfer functions for the in-ground-effects case are obtained 
by including the terms , C , and C , as presented in the previous section, 

"(j) >(() 

where bank angle, (}), is considered proportional to the difference in height 
between the right and left v/ing. A change in height produces a change in later- 
al moments and accelerations. This, in turn, would be reflected in the lateral 
transfer functions, as indicated in table 5. The root loci of the two aircraft 
control system configurations with outer loop closure (K^ increasing) show 

similar critical gain and frequency characteristics. The equivalent character- 
istic poles of the lower modes (lu^ = 1.25) with the attitude loop open indicate 

a lower damping ratio for the aircraft in ground effect, as shown by the dashed 
slope of figure 9. The same trend is also evident in comparing the two charac- 
teristic equations of the basic aircraft in and out of ground effect (table 5). 
However, this mode (u)^ = 1.25) for both configurations is heavily damped, and 

any difference thus shown would be unnoticeable by the pilot and not expected 
to have contributed to the PIO problem. 

Nonlinear calculations. - During the PIO, the flight control system was 
extensively rate and position saturated, as evident in the time history shown 
in figure 3. Therefore, a realistic analysis requires consideration of these 
nonlinearities. The discussion that follows in this section utilizes a non- 
linear computation of system gain as a function of frequency. This is done by 
making various amplitude and phase comparisons of four conditions: the initial 
system on the prototype aircraft and the MFCS both in and out of ground effects. 
Two pilot inputs are shown for comparison in subsequent figures; F = 4.448 

newtons (1 pound) and F^ = 35.68 newtons (8 pounds). The lateral stick force 

per roll rate varies according to the feel gradient shown in figures 5 and 8. 

The method is simply to assign a value to the error vector (amplitude and 
phase). Then an input/output vector can be determined through each element, 
keeping track of amplitude and phase algebraically until an open-loop vector 
is obtained. This makes any vector or vector ratio a function of input and 
frequency. This procedure is repeated on the digital computer a sufficient 
number of times to saturate the roll command for the frequency range of inter- 
est. For low inputs within the rate and position limits of the system, a 
linear calculation is used as a computational check. 
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Appendix B presents details of the nonlinear computational method. The 
aircraft transfer functions are presented in table 5. 

Amplitude and phase plots for the initial aircraft control system, which 
existed on the prototype aircraft at the time of the first flight, are pre- 
sented in figures 10(a) and 10(b), respectively. The attenuation is normal 
with frequency. Higher lateral forces reflect more system gain, due primarily 
to the nonlinear feel gradient. However, high forces add more phase lag due 
to rate and position saturation. Above 2.6 radians per second a force input 
of 35.58 newtons (8 pounds) will rate saturate the system throughout the higher 
frequency range. Near the PIO frequency (lu = 4.0 rad/sec), the calculations 
indicate a phase change of -180 degrees. 

The initial system in ground effects (figs. 11(a) and 11(b)) in general 
shows reduction in system gain at the lower frequencies, with little differ- 
ence evident at the high frequencies. The crossover frequency and critical 
gain margin are about the same as those shown for the out-of-ground effect 
configuration. 

Figures 12(a) and 12(b) present an amplitude and phase variation, with the 
frequency of the MFCS configuration in ground effect over the same frequency 
range as the initial system on the prototype aircraft. The system gain has been 
reduced sufficiently so that the system behaves in a linear fashion for the two 
lateral force inputs shown. This is apparent if it is noted that the gain dif- 
ference between the two force levels is the same and that the phase angle vari- 
ation is the same throughout the frequency range. 

In order to sustain an oscillation in a closed-loop system, the total gain 
must be unity at a phase shift of -180 degrees (minimum conditions, according 
to the Nyquist criteria). If it is assumed that the pilot acts as pure gain, 
with no phase changes, the pilot gain required to cause a PIO would be the gain 
margin at a -180 degree phase shift. Because of the nonlinearity, the critical 
pilot gain would also be a function of the lateral stick force input, F^. 

Figures 13 and 14 present PIO frequency and critical pilot gain variation as a 
function of lateral stick force for various control configurations. Low force 
input approaches the linear values of crossover frequency and gain as predicted 
by the root locus (fig. 9). 

Ground effects, as noted previously (fig. 9), have little effect on the 
critical pilot gain. For reasonable values of input, the critical pilot gain 
with the MFCS has been increased by a factor of 4, or at least 10 decibels. 

This is indicated by the level of the two curves (F /(j) ) as a function of F 

3 0 cl 

shown in figure 14. Also, as is evident by the constant value of PIO frequency, 
the MFCS configuration behaves in a linear fashion throughout the force level 
presented. 

A detailed investigation of the initial control system using the nonlinear 
calculations indicates that rate saturation exists throughout the higher fre- 
quency range (from 2.6 to at least 8 rad/sec) at the highest command response. 
Total flaperon surface position difference versus input frequency is presented 
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in figure 15 for the initial control system in ground effect. Rate limiting 
(Aio =112 deg/sec) is shown by the hyperbolic boundary. The command response 
to the flaperon is position limited at the lower frequencies and decreases 
according to the first order lag of the actuator. As noted, above 2.6 radians 
per second, flaperon surface is limited in travel by the maximum rate. 

A comparison of flaperon surface limitations is presented in figure 16 
for the MFCS configuration for three levels of force input. A force input of 
35,6 newtons (8 pounds) causes no surface saturation. A force input of 53.4 
newtons (12 pounds) will be marginally rate limited at high frequencies. Above 
71,2 newtons (16 pounds), maximum roll rate is commanded and the system will 
rate saturate above 3.5 radians per second; this is below the flaperon surface 
1 imit. 


CONCLUSIONS 


A lateral pilot-induced oscillation (PIO) experienced on the first flight 
of the YF-16 aircraft was analyzed by linear and nonlinear techniques. In 
general, the analysis indicates: 

(1) The nonlinear technique can be used to analyze this type of control 
problem. For small inputs, where saturation does not occur, the linear and 
nonlinear computations agree. 

(2) The pilot system phase relationship (bank angle/pilot input) during 
the PIO was at least -180 degrees. The PIO frequency for out-of-phase response 
was 3.2 radians per second. Under high pilot gain conditions, the initial 
flight control system response and control actions could result in a PIO as 
experienced on the first flight. 

(3) The rate command system essentially behaves the same dynamically 
regardless of the proximity of the ground plane. 

(4) Control system modifications (MFCS) for the takeoff and landing 
phases have increased the pilot gain required to PIO the aircraft by about 10 
decibels and also made the mechanical system more difficult to position and 
rate saturate. 


Dryden Flight Research Center 

National Aeronautics and Space Administration 
Edwards, Calif., August 24, 1979 
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APPENDIX A - GENERAL DYNAMICS FLIGHT REPORT 


The first flight of the YF-16 airplane was conducted at Edwards, Califor- 
nia by General Dynamics Engineering Test Pilot P. F. Oestricher. The follow- 
ing material is the flight test report in its entirety as submitted by the 
pilot, and represents a rather complete description regarding the pilot's 
impression of the PIO experienced during the unintentional lift-off. 

"The purpose of these series of tests was to perform a limited functional 
check of various systems (including the instrumentation system and test control 
at Bldg. 3940 and the trailer) and to determine the taxi characteristics at 
various speeds. 

"The test configuration was that of the basic airplane with an AIM-9 
missile mounted on each wingtip. The airplane was fully fueled at the start 
of the tests and was ready in all respects. 

"Taxiing at normal speeds was evaluated while moving the airplane to the 
"last chance" check area for runway 22. Periodic application of brakes was 
requiredto prevent an excessive speed buildup. The braking effort expended 
by the pilot (product of pedal force and duration of pedal displacement) was 
perhaps 30 percent to 50 percent more than required in the case of a fully 
fueled, clean configured RF-4C. Nose wheel steering was used throughout the 
runand proved to be precise and easily controlled. Following a check by the 
mobile crew, the airplane was positioned on runway 22 for an idle power taxi 
run without brake restraint. A taxi speed of around 30 knots was noted during 
this test. After a period of straight ahead taxiing, several S- turns were 
made with no difficulty. The airplane was stopped after traversing about 
5,000 feet. 

"Following an inspection by the mobile crew, the airplane accelerated 
toward a target speed of 80 knots. It is believed that an overshoot of about 
10 knots occurred on this run. The nose wheel steering appeared to be overly 
sensitive at speeds of 50 knots or higher and was accordingly disengaged. 
Directional control by rudder was very satisfactory after the NWS disengagement. 
The airplane was stopped using moderate brake pedal force after traveling about 
5,000 feet. It was then towed back to the "last chance" area for runway 22 for 
brake cooling. 

"The brakes were checked and found to have cooled sufficiently to resume 
the taxi tests. A normal start was accomplished as were the pre-takeoff check 
list items. The IIRS was aligned and checked for proper operation. The air- 
plane was positioned on runway 22 for the planned 135 knot high speed taxi run. 
The brakes were held and the power level slowly advanced to determine the RPM 
at which wheel slide would occur. This was determined to be about 87% RPM. 

The gross weight at this time was about 21,200 lbs. The corresponding C.G. was 
34.3% M.A.C. The engine was kept at idle RPM until the runway winds (as 
reported by the tower) dropped below the 12 knot maximum agreed to for the taxi 
run. Upon tower clearance for the run, the brakes were released and intermedi- 
ate power selected for a period of about six seconds, after which a substantial 


16 



power reduction was made. Nose wheel steering was disengaged at an estimated 
50 knots. At about 130 knots (but apparently with the airplane still acceler- 
ating somewhat) the airplane was rotated to about 10 degrees and small lateral 
stick inputs were made in an attempt to get a feel for control response. No 
response was noted by the pilot (doubtless because the main gear was still 
restraining the airplane from rolling) and the angle of attack was intentionally 
increased a small amount. The airplane had continued to accelerate during this 
time but the pilot was unaware of the fact. Immediately upon rotating the 
second time, the airplane lifted off with the left wing dropping rather rapidly. 
Right roll command was applied and the airplane was immediately involved in a 
fairly high frequency lateral PIO (10 cycles in 14.3 seconds). Eventually, the 
roll oscillation was stopped but not before lightly touching the rolleron wheel 
on the lower outboard fin of the left AIM-9 to the runway, striking the right 
horizontal tail tip (at the trailing edge) on the runway, bouncing off the main 
landing gear several times in a nose-high and generally symmetrical manner, and 
developing a substantial heading deviation from the runway axis. The latter 
factor promptedthe decision to fly out of the situation as it was felt that it 
would be impossible to steer the airplane so as to remain on the runway, even 
if the nose wheel could be quickly brought down to the surface. Intermediate 
power was applied for a short period of time after which a fairly low thrust 
level was held. The airplane was allowed to slowly climb away in a shallow left 
turn, with a minimum of pilot control inputs being made. A downwind leg to 
runway 22 was established at about 600 feet AGL at 175 KIAS. The ADC caution 
light was noted to be on at this time. No attempt was made to turn the light 
out by resetting. A wide pattern was flown to a long, decelerating final 
approach with 12 degrees being established just prior to touchdown. A very 
slight (low amplitude and frequency) lateral motion was noted prior to touch- 
down. The ground effect was quite pronounced and the engine was brought to 
idle while still airborne. Aft stick force was relaxed after touchdown and 
the nose wheel fell gently to the runway, at which time the speed brakes were 
commanded open. It should be noted that the pitch trim was still in the 
neutral position at landing since no pilot trim had been applied during the 
flight, _ Moderate braking was applied until the airplane was stopped. Follow- 
ing an inspection by the mobile crew, the engine was shut down and the airplane 
was towed to the hangar. 

"The tactics attempted during the PIO are evident from watching the excel- 
lent movie films available. Briefly the attempt was to: 

1. Keep the wingtips off the runway and stop the roll oscillation with 
the wings level . 

2. Recover the nose high attitude when the lateral control problem had 
been solved. 

3. Control altitude and vertical velocity with thrust. It is believed 
that this particular attempt was relatively successful. 

"No sideslip was noted by the pilot at any time despite the violent nature 
of the oscillation and the full lateral commands being applied. The roll con- 
trol problem appeared to be the most serious by far, and accounted for most of 
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the pilot's attention at the time. Once away from the ground, and the need to 
keep roll angle within tight bounds, the pilot was able to relax, with the 
results which are evident in the movie film. The pattern and landing were 
understandably somewhat conservatire [sic] although a small rudder doublet was 
performed during the final portion of the approach in an attempt to assess 
directional control sensitivity. No dihedral [sic] effect was noted and the 
airplane felt somewhat sensitive as compared to other tactical airplanes. 

"Takeoff and landing gross weight/C. G. combinations were 21,100 lbs./34.3% 
M.A.C. and 20,300 lbs./35.0 % M.A.C. , respectively." 
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APPENDIX B - BASIC NONLINEAR ANALYSIS AND COMPUTATIONAL PROCEDURE 


Nonlinear Analysis 

The functional block diagram shown in figure 17 is used for total lateral 
control system nonlinear definition and analysis. Each and all forward, feed- 
back and common paths are as indicated or represented equivalently. The initial 
task is to further separate the common paths through the rudder loop so that 
the overall mathematical description of the open-loop attitude system will con- 
tain only one summation junction. This finalizeddescription will then enable 
an explicit solution for the input by adding the roll rate feedback (K p) to 

the error signal commanded ((})g). 

From the figure 17, the following set of equations are evident. 


= (Vy * “i'"' ■ V * ♦e) 


where 




and 


<5p - Gp <|)g 
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The nonlinearities are appropriately included in subsequent numerical 
computations. Substituting A , r, and p into the 6 equation and rearranging 
gives: ^ ^ 


r _ 



a 


Vf ^ ^ 


—HU) 


H Gr ^ 

P F 6 r 


( 1 ) 


where 


H(s) = 1 - HyG, - H,G, ^ . HpG, §- 


(2) 


The total roll rate commanded per error signal (|)g is: 


<l>e ^ If' 4>e 


(3) 


or 


1 m) — 6" r 

<}>e ' r (})g 


(4) 


Substituting H(s) and into the above equation, and then multiplying through, 

^e 


canceling and regrouping the aircraft transfer functions with respect to roll 
response yields: 


£ 


{h * °F - Vr ij • 







•< 

'I + G ^ 

•j •'V 

‘ H, (< 

\ A 

^ ^ ‘^F 6 

a 

(' 

r , p r 

’H 63 6, 


m 


(5) 


Since no nonlinearities are present in the rudder closed- loop feedback 
function H(s), this function in subsequent computations can be conveniently 
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moved downstream to the roll rate and roll attitude loops as and 

respectively. The aircraft open-loop transfer functions are listed 

in table 5. The actuator transfer functions Gp, and G^ (including their 

effective gains) and the feedback functions H , and H are as follows: 

V r ’ D 


r - 

■ s + 20 

o - 20 

" s + 20 


Hy = -296 

u _ 1.33 (3s + 15)s 
^r (s +15) (s + 1) 


_ 11.25 
r s + 20 


_ .345 (3s + 15)s 
p (s + 15) (s + 1) 


In order to list and define the flaperon rate and position limit boundary in 
computing the total roll response, equation (5) is separated so as to reflect 
the output of each individual surface. 


L 

({) 

®(ai leron) 




r 6^ 





( 6 ) 


^e 

(flaperon) 



H G ^ - H G ^ 
y r 6^ r r 0 ^ 


(7) 


(j) 

®(rudder) 


r 6 ^ 




+ H. 




( 8 ) 


The roll contribution is small due to and r in the above rudder equation 

(eq. (8)). However, for the sake of completeness, all terms are included in 
subsequent computations. 
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Computational Procedure 


The flow diagram shown in figure 18 represents the YF-16 aircraft roll 
control system, aerodynamics, and equivalent rudder interconnect and feedback 
loops. 

In the diagram, T1 represents the pilot as a pure gain constant. Thus, 
the commanded input is proportional to the error in bank angle X 2 

(X 2 “ ‘I’r " output of any quantity X^, X^, and so forth de- 

pends on X 2 and the intermediate transfer function T2 and T3. T2 is the 
transfer function for each of the feel systems shown in figures 5 and 7. 

For both system configurations, T3 is the forward loop model of the 

dual lag system. 

As previously shown, it was convenient to rearrange the rudder inter- 
connect loops such that the functions T9, T13, T16, T17, and T18 gave an 
equivalent representation for the separated loops. The aerodynamic transfer 
functions for both in and out of ground effects are listed in table 5. In 
order to express the subsequent rate limit Til as the total difference between 
the two surfaces, the gain constant of T13 is equal to half the listed value. 
This representation can be used for both linear and nonlinear calculations. 
However, for nonlinear calculations the input/output ratio, both amplitude 
and phase, will be a function of the input. 

To obtain a frequency response, various values are assumed and assigned 
to Xg (A, lu), a vector that is a function of both frequency and amplitude. 

Then the vector signal is traced for each calculation. The phase and ampli- 
tude of each successive vector are noted through each block and summation 
junction until an input/output vector relationship is determined. 

The logic, vector algebra, and transfer functions shown in table 6 for 
various configurations were programmed on a digital computer. In the equations, 
s = ju), where w is a variable, and Xg is a variable. 
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TABLE 1. - PHYSICAL GEOMETRY AND AIRCRAFT CONFIGURATION 

s = 26.0129 (280 ft^) 

c = 3.334 m (10.937 ft) 

b = 8.891 m (29.17 ft) 

w = 9616.157 kg (21,200 lb) 

c.g. = 34.3 percent 

I^ = 12472.57 kg-m^ (92,200 slug-ft^) 

I^ = 78631.43 kg-m^ (58,000 s1ug-ft^) 

1^2 = 203.36 kg-m^ (150 s1ug-ft^) 


TABLE 2. - YF-16 FLIGHT CONTROL LIMITS 
(a) Control surfaces 



Authority limits 

Rate limits 

Horizontal tail 

±25 degrees 

60 deg/sec 

Flaperon 

±20 degrees 

56 deg/sec 

Rudder 

±30 degrees 

(b) Command servos (in series) 
Second order characteristics 

120 deg/sec 

Horizontal tail 

C = .7 a. = 52 

175 deg/sec 

Flaperon 

C = .7 & = 52 

150 deg/sec 

Rudder 

C = .7 & = 52 

210 deg/sec 


(c) Secondary surfaces 

Authority limits 


Leading edge flaps 0 30 deg 

Trailing edge flaps 0 20 deg 

Speed brakes 0 60 deg 



TABLE 3. - FLIGHT CONDITIONS 

a = 15° 

M = .20 

h = 694.94 m (2280 ft) 

r 

q = 2614.248 N/m^ (54.60 Ib/ft^) 
LEF = 25° 

T.E. flaperon (roll control) < ± 20° 
V = 67.57 m/sec (221.7 ft/sic) 

6o = 15° 


TABLE 4. - REVIEW OF RELEVANT PIO CONDITIONS 


(a) First flight of YF-16 

1. Frequency =4.2 rad/sec, .668 hz 

2. 12 complete cycles 

3. Command input « 160 percent of saturation 

4. Total lag between pilot force input and bank angle -output > 180° 

5. Flaperon rate limit = 112 deg/sec (= 2x surface rate) 

6. Horizontal tail rate ■ less than system rate saturation 


(b) Conditions common to PIO problems 

1. -180° pilot and system relationship 

2. Reasonable frequency range 

3. Rate limiting 


(c) Differences from other aircraft with PIO problems 

1. Minimum displacement force stick 

2. Ground effect (takeoff) 

3. Fly-by-wire (no mechanical overlapping loop) - consequently no "stick talk" 

4. High reponse roll command (design requirement) 
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TABLE 5. - AERODYNAMICS 


(a) Body axis stability coefficients 




-.00508 


-.001 

\ 


.0035 


-.00027 

s 

* 

-.019 


0. 



-.205 


-.0023 

"P 


-.08 

s* 

-.00024 



.21 


0. 



.198 


. 00047 



-.54 


-.00173 


— 

.9 


. 00278 
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TABLE 5 - Continued 

(b) Transfer functions - no ground effect 


£ = “2.77(5 - .038)(s + .2ZS ± il.438)^ 


£ - -6.35(s - .038)(s + .217 ± 11. 308V 

a 


1.29(s - .Q39)(s - 2.63)(s 2.50) ^ 

A 


126(s + .528)(s - .344 ± .j3.52V 
A 


122(s + .524)(s - .976 ± .j4.88V 
A 


757(s + ,559)(s + .122 ± il.69V 
A 


= -.049(s - 1.286)(s + .384)(s . 357 ± il.90^ ^ 

A 


= -.047(s - 3.233)(s + .328)(s .730 ± 11.80) ^ 

A 


_ -.180(s - .108)(s + .567)(s + .136 ± i2.48)^ 


where 


A = (s + .719)(s + .093)(s + .114 ± j2.27)‘ 
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where 


TABLE 5. - Concluded 

(c) Additional terms to include ground effects 

C. = -.0260 
^«l> 


= .00635 


n. 


' 4 > 


= 1.376 

y<!> 

(d) Transfer functions including ground effects 
E _ -2.77(s - .409Ks .227 ± 11.428) ^ 


E = -6.35(s - .0404)(s •*• .218 ± il.301) ^ 


E = 1.29(s - .043)(s - 2.59KS 2.47V 


r _ -.126(s .478Ks - .319 ± 13.82) ^ 


I = -.122(s .485)(s - .957 ± i5.23) ^ 


r = -.757(s .436)(s •*• .184 ± 12.00) ^ 


= -.049(s - .297 t ,11.008)^(5 ^ .214 ± 12.33) ^ 

A 


= -.047(s - 1.11 ± .1.736)^(s + .390 ± i2.31) ^ 

A . 




•.180(s ■>••.210 ± i.770)^(s + .155 ± 12.53) ^ 


A = (s + .373 ± j.366)^(s + .144 ± j2.48)^ 
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TABLE 6. - TRANSFER FUNCTION DEFINITIONS 
(a) YF-16 prototype aircraft 


Transfer 

function 

T1 

T2 

T3 

T4 

T5 

T6 

T7 

T8 

T9 

TIO 

Til 

T12 

T13 

T14 

T15 

T16 

T17 

TIB 


Prototype ' 
Change 


Transform 

(fig. 18 ) 

X3 = X2*l 

X4 = .008726(X3)^ + .06980803(X3) 

If )X4| > 3.5, X4 = (3.5/|X4| )*X4 
X5 = (10/(s + 10))*X4 
X7 = X6*.2 
X8 = .25*X7 

X9 = ((20*s)/(s + 20))*X8 
XIO = X9*l 

If iXlOl > .934, XIO = (.934/|X10| )*X10 
xn = XlO/s 



X13 = ((20*s)/(s + 20))*X7 
X14 = X13*2 

If |X14| > 1.868, X14 = (1.868/ iX14( )*X14 
X15 = X14/S 

If |X15| > .698, X15 = (.698/|X15| )*X15 

X16 = ^yX15 

X17 = .563*X7 

X18 = (20/(s + 20))*X17 




*X18 


X22 = 




/r 




r r 


r 








(b) Modified flight control system 


:ontrol 


T2 


system 


T4 


MFCS 

to X4 = ,0007487(X3)^ + .01745(X3) 

|X4| <2.91 
X7 = X6*.12 
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Figure 2. Simplified roll and yaw axis flight control system. 
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Figure 3. YF-16 pilot induced oscillation ex 

perienced on the first flight. 
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Figure 5. Initial roll command gradient on the prototype 
aircraft. 
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Figure 7. Diagram description of 
mechanical system. 
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Figure 9. Root locus of altitude loop for the 
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Figure 10. Concluded. 
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Figure 11. Concluded 
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Figure 14. Comparison of the PIO frequency and minimum 
pilot gain in ground effect of the initial" control con- 
figuration to the MFCS. 
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Figure 18. YF-16 mathematical model 
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